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^ ■ INTRODUCTION 

■ , Random Matrix Theory (RMT) has been a useful laboratory for simulating hamil- 

Q I tonians of statistical systems with random matrix elements^, due possibly to impurities 

>-^ ■ scattered around the material. An unphysical yet simplest choice is to assume each of 

matrix elements of the hamiltonian M to be independently derived from the gaussian 
^^ . distribution, 

^: Z= f dMe-^*=-^^', (1) 

JNxN matrices 

to which is associated Wigner's single-banded, semi-circle spectrum. This macroscopic 
spectrum is, however, by no means realistic; physical spectra are far more complicated. 
Namely they may even have multi-band structures, which in the side of RMT corre- 
spond to non-gaussian distributions. Thus, it is only universal quantities insensitive 
to a chosen matrix measure that may be justifiably extracted from RMT for physical 
systems. 

Various quantities concerning microscopic spectral correlation in the bulk of the 
spectral band have long been believed universal, as identical results are derived from 
non-gaussian ensembles corresponding to classical orthogonal polynomials other than 
Hermite^. Here the term 'microscopic' refers to that the correlation is measured in the 
unit of the mean level spacing, which is of order 0{1/N) for ensembles of the type (|l]). 
This conjectured universality of the sine kernel 

sin7r(A-A0 
^^^^'^)= 7r(A-V) ' (2) 
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which comprises all spectral correlators 



p,(Ai,---,Ap) = (n4*r5(A„-M))= det Ks{\a,Xb), (3) 



, -, N / l<a,b<p' 

\a=l I ~ 



is proved^' ^ for unitary invariant ensembles with generic single-trace potentials e"^*''^*-*^/ 
RMT has recently extended its range of applicability toward QCD^. There the 
Dirac operator ip = 7^(z 9^ + 74^), having gauge fields as random elements, is regarded 
as a random hamiltonian. Schematically, the Euclidean QCD partition function^ 

Zqod = /|<i^,«,*,|e-?J''"''-/'''-^'"''*' 

is transformed by the change of integration variables 

A,^iJI)a^M={1^ ^M (5) 



into a RMT 



^RMT 

/Afx TV matrices ^T \ W Ifl'f 



f df,{W)Y[detl "^ ^ ]. (6) 

J Afx TV matrices r \ V^ ^f / 



Here N stands for the size of the Dirac operator i.e. the volume of the spacetime, and 
d^{W) a measure invariant under 

W ^UWV^ {U,V e V{N) for W complex). (7) 

Novel features of this RMT in addition to the conventional model (|I|) are: 

• Presence of the fermion determinant. For small rrif, it expels Dirac eigenvalues 
of order 0(771 f) from the origin. 

• Chiral structure of the random matrix. Since {M, 75} = each eigenvalue A is 
accompanied by its mirror image —A. The Coulomb repulsion between these pairs 
prevents them from populating around the origin. 

Among the quantities calculable from this type of RMTs, of particular interest are the 
distribution and correlation of soft eigenvalues of the Dirac operator as their accumu- 
lation is responsible for the spontaneous breaking of the chiral symmetry via^ 

s = m)\ = ^. (8) 

Due to the above two features, correlations of soft eigenvalues are expected to deviate 
from the universal sine law (0). Verbaarschot and his collaborators^ have calculated 
correlations of eigenvalues of order A ~ 0{1/N) from the gaussian ensemble 

dfi{W) = dW e-'''''''^''^ , (9) 

appealing to the conjecture of universality. These analytic predictions have been com- 
pared first to the numerical data from the simulation of the instanton liquid modeF, 
and more recently to that of the lattice gauge theory*^' ^. The agreements are impressive. 



'For simplicity, we consider only the topologically trivial sector. 



In this article I shall review the proof of universality of the microscopic spectral 
correlations of the model @ for rrif = and with single-trace potential measures^, in 
order to justify partly the use of gaussian ensemble (^) and to corroborate the above 
mentioned agreements. Rather than to follow closely the original proof ^^' ^"^5 I shall 
disguise it with the conventional approach of Q, P operators^^ and with the method 
employed in a recently proposed alternative proof ^^' ^^, in a hope that borrowing notions 
from quantum mechanics may add the proof some pedagogical flavor. At the end of 
the article I shall address a problem associated with multi-criticality. 

ORTHOGONAL POLYNOMIAL METHOD 

We start from recalling basic technology of random matrices: the orthogonal poly- 
nomial method. We consider the chiral unitary ensemble (xUE) 

^xUE = f rfMe-^-^(*^^) IdetMr, M = f ° ^M (10) 

capturing the global symmetries oi Nc > 3 QCD4 with Nj = a massless fundamental 
fermions, as well as the unitary ensemble (UE) without chiral structure 



>UE 



/ dMe-^*^^(*^') IdetMr (11) 

J NxN hcrmite 

modeling QCDs^*^. We allow non-integer a > —1, for they can be treated on the same 



footing as integer cases. Both ensembles allow eigenvalue representations^^ 

/■oo ^ , ^ ^, 2 

Z^UE = / Yl (dXiXte-^^^^^n dety-^ , (A^ : eigenvalues of VT+l^) (12 
JO ~ri ^ ^ ^3 



coo " 

^-00 .^-^ 



det A, 



i-l 



(Aj : eigenvalues of M) (13) 



Inside the integrals, det A* ^ may be replaced (up to an irrelevant constant) with poly- 
nomials detPj„i(Aj) orthonormal with respect to the measure dX |A|"e~''^^'''*'^: 

5^m= /rfA|Are-^^Wp„(A)P„(A)= /dA^„(A) ^^(A). (14) 



Here we have introduced one-particle wave functions of a free fermion (at the n-th 
state) 

V.„(A) = |Ar/V^^/^P„(A). (15) 

The fermionic nature comes from the Vandermonde determinant. The wave functions 
at first N levels comprise the Dirac sea, i.e. the A^-particle ground-state wave function 

vi>^({A})= det V^„_i(AJ. (16) 

l<n,m<N 

Using this ground-state wave function, the vacuum expectation value of unitary invari- 
ant observables is written as 

{0)= jd''X^>N{{X})0^>M{{X]). (17) 



"i- Microscopic universality for the gaussian potential plus certain terms which breaks the invariance (uf) 
is proved in refs.lO, 11. See also ref.l2 for a perturbative treatment of a quartic potential. 



Namely, spectral correlators are related to tpn by 

/ p 1 



\ Wick's theorem 
p^(Ai,---,Ap) = (n^tr<5(A,-M)) ^ det K^{Xa,Xk) 

\a=l N I ^^»'^^P 



where Kn is the projector to the Dirac sea: 



N-l 



Kn{\\')= E^"W^"(^') 



n=0 



i<N(\)i<N-l{\') - 'll>N-l(\)i<N(\') 



Use is made of Christoffel-Dalboux formula in the last line, and q^q is a constant defined 
immediately below. 

We note that the orthogonality relation (JT^) for the xUE case (/q°° dX) can absorbed 
into the UE case (/^ dX) by the change of variables A — > V^ accompanied by the 
redefinition a ^ ^ (or a ^ ^) and Psnl^A) ^ Pn(A) (or yAP2n+i(yA) ^ 
P„(A)). Thus in the following we need only to consider the UE case, where orthogonal 
polynomials have definite parities. 

Now let us for a while concentrate on a = 0, where the logarithmic component of 
the potential is absent. In this case one may represent the Heisenberg algebra [^, A] = 1 
by matrices Q and P acting on the Hilbert space of orthonormal wave functions^^, 

A V'n = gn+lV'n+1 + ^n^n-l = JI Qnm ^m, (20) 



d\ 



Ipn 



N 



V ipn + (linear comb, of iJn-i, i'n-s, ■■■) = Y1 P"™ V^™- (21) 



(Anti-)self adjointness of multiplication and differentiation operators inherit to the 
(anti-) symmetry of the representing matrices. 



d\ {\iJn)i^n 

dX I -TTlpn j -ipm 



dX i)n (A iJrn) 
- I dXtpn (tV^™ 



P« 



_P 

-■- m.n 



(22) 
(23) 



Although eq. (^TD and the antisymmetry (|2^) are enough to determine the P matrix as 

Pr,™. = T— V'CQ)™™, (n.^m), (24) 






we shall not use this fact except for fixing a integration constant in the sequel. 

At this stage we consider the situation when the coefficient qn become single- valued 
in the limit n, m ~ A^ ^ 1, 

Qn ~ qn±i ~ g = 0(1), (25) 

a crucial assumption which is violated in matrix ensembles with multi-band spectra^". 
Then the matrix elements of Q, P for ra, m = A^ -|- 0(1) take the forms: 



/ 



Q. 



\ 



q 
q q 
q q 



I 



N 



■ J 



••. ••. ••. ••. \ 

■•.■•. -P2 

■•• -pi -P2 

'■■Pi -pi 

P2 Pi -pi 

P2 Pi '■■ / 



(26) 



A matrix of the above form of P can always be written as 



^ fc>0 ^ 

where the A matrix stands for the antisymmetrizer 

/ . . \ 



A, 



(27) 



V 






-1 




1 





— 




1 




1 



:28i 



■• / 



When a 7^ 0, ^ can not be represented on the space of orthogonal polynomials as 
it involves a term proportional to 1/A. However one may verify^^ that the operator 
1\ ~~ 2 A " ^^ indeed represented by the same P matrix as in (P7|). Here n is the 
number operator nipn = nipn, and the existence of the signature (— )" preserves the 
antisymmetry of P, 



V'2m+l ((-)"V^2n) = - ((-)"^2m+l) i>2n- 



(29) 



In the next section eqs.(^) and (27) are utilized to prove the microscopic universality. 



PROOF OF UNIVERSALITY 



Here we re-exhibit the asymptotic forms of recursion relations (p6D, (pT]) satisfied 
by the orthonormal wave functions (a = 0): 



Xlpn = q {ipn+l + i'n-l) 



^ / 



NA{X) 



{-^Pn+l + ^n-l) • 



(30) 
(31) 



They are combined to yield unharmonic analogues of lowering and raising operators of 
a harmonic oscillator at highly excited levels n ^ 1: 



^NA{X) dX 2q 
1 d X 



+ 7r]i'n = Vn-l, 



NA{X)dX 2q 

'W'n 



+ TT ] i'n = i'' 



n+1- 



(32) 
(33) 



The requirement that two processes ^n\t"^^J''ipn niust commute amounts to the sup- 
pression of the commutator of the raising and. lowering operators in the large-n and A^ 
limit. 



A 



d 



Oi 



<1. 



2g ' NA{X)dX\ ' N* 

Thus we are lead to a second-order differential equation for ipn, 

2 / , \ 21 



d 



NA{X) dX 



1 - 






-ipn 



0. 



(34) 



(35) 



As mentioned in the end of the last section, the only modification for a 7^ case 
is the replacement 



A/ A. 

Z\^" ^ [dX 

The resulting equation takes the form 



a 



2X 



i^n- 



1 1 


f d (-)"a^ 
[dX ' 2A J 


1 
A{X) 


f d 
\dX 


i-ra^ 


+ 1- 


(X 


N^ A{X) 


2X ) 


\lq 



4^n 



Let us first take the macroscopic (large- A^) limit: 

N ^ 00, X : fixed 
of (|37|). Since 1/A^ plays the role of h in this limit and thus eq. 

1 1 d^ 



is reduced to 



m A{xy dX' 



+ 1 



A 

2g. 



^n = 0, 



(36) 



(37) 



(3^ 



(39) 



we obtain the WKB solution'^ |: 



ipniX) occos A^ / dXA{X) 

V 



'\ 



A 

2g. 



+ 



nn 



(40) 



Substituting this solution to (0), we obtain the kernel in the large- A^ limit: 



K{X,X')= lim Kn{X,X') oc 



sml^NJ^,dXAiX)^l-{^^) 
A -A' 



\2 



p(A) = i^(A,A) = const.A(A) 



\ 



[^1. 



(41) 
(42) 



The constant is not fixed within this approach, though it can fixed to be l/vr if we use 
the explicit form of the P operator (P^). The meanings of q and A{X) become clear 
at this stage: ±2g stand for the edges of the spectrum and A{X) the deviation of the 
spectral envelope from Wigner's semi-circle. 

Now we proceed to take the microscopic limit: 



A^ 



00, A ^ 0, z = NX : fixed. 



(43) 



We have already assumed the single- valuedness of the recursion coefficients g„, which 
is known to be the sufficient condition for A{0) = vrp(O) > 0. (In other words, we have 
assumed the chiral symmetry breaking as an input.) Then eq.(^7|) is reduced to 



d_ (-Ya 



A 
dz 



"a 



.4(0)2 \dz ' 2z J \dz 2z 

Its solution which is regular at 2; = is a Bessel function: 

/ e\ 

ipn OC \G Jail {A{0)z) . in = 

2 V o 



+ 1 



^n = 0. 



even 
dd 



(44) 



(45) 



§ This process was initiated in the context of stochastic quantization of matrix niodels^^, ahhough 
their Fokker-Planck equation differs from eq.(P^. 1 thank C. F. Kristjansen for remarks on this point. 
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Within either even or odd sector of the wave functions, the n- dependence enters only 
through An{0). Thus we have 

i^n - ^n-2 ~ ^^V-n OC ^i (j^ {A{0)z) - J^ (A(O)z)) . (^U = 'J^) (46) 

Substituting the above solutions again into ([19D , we obtain universal forms of the 
microscopic kernels (called the Bessel kernels^^) 



X-UE: A' 



(.. .') = ,^^o^^^-J^.d2Mo>-)JfMoy)-{^»^l ,,3) 



2;^ — z' 



2 z — z' 

and microscopic spectral densities Ps{z) = Ks{z,z), 

XUE : psiz) = inp{0)f \z\ (j^ - Ja+iJa-i) (2vrp(0)^), (50) 

UE : ps{z) = l^^r^\ z fjl+i + JLi - J^+^J^^ - J^^Jc^ (7rp(0)z) . (51) 

In the above, the integration constants are fixed by requiring \m\z^QoPs{z) = p(0). 
When measured in the unit of mean level spacing l/p(0), Ks{z, z') dz contains no free 
parameter. After this rescaling, the Bessel kernels (^8]) and (^91) approach the sine 
kernel (0) in the limit z, z' -^ oo, z — z' = 0(1) as the repulsive effects due to the 
fermion determinant and the chiral structure become negligible. 

DISCUSSIONS 

In this article I have presented a proof of universality of microscopic correlations 
for RMT modeling QCD3 4 coupled to massless quarks, within a single-trace potential 
class. The spectral correlation function is shown to be insensitive to the choice of 
matrix measures as long as the macroscopic spectrum is supported on a single interval. 

Universality of the Bessel kernel is valid even beyond large- A^ matrix models. 
Namely, the zero-dimensional reduction of the SU(A^/) a-model of pions also yields 
sum rules^^ which are identical to the moments of (0). Together with the numerical 
agreements mentioned in the introduction, it is very convincing that this wide range of 
universality encompasses QCD4. 

Recently it also became possible to incorpolate quarks whose masses are within the 
same microscopic range as the Dirac eigenvalues^^' ^^. We expect that the exactness of 
those universal correlations be verified by numerical simulations in a foreseeable future. 

On the other hand, one may address a question: whether RMT can describe chiral 
symmetric phases of QCD such as Higgs phase or the large Nj case. In these cases, 
accumulation of the soft eigenvalues is not strong enough to form a chiral condensate, 
although the theory is still interacting. There is no satisfactory answer to it so far, yet 
one may check the existence of universality within RMT. If the microscopic universality 
in the vicinity of the origin is broken for random matrix ensembles with p(0) oc | (ipip) \ = 



0, we may not expect to extract informations from RMT. In view of the above, we are 
lead to analyze a RMT with a double-well potential such as 

Z= f dMe-^*^(-^'+^^'), g = - (52) 

JNxN matrices 4 

whose coupling is tuned so as the macroscopic spectrum p(A) to have a double zero at 
A = ^^. According to ( ^2]) it corresponds to A{X) of the form 



A(A)~A2 + 0(A^). (53) 

Then eq.(|37D suggests a possible scaling limit: 

iV -> oo, A -^ 0, z = N^/^X : fixed. (54) 

However, the whole procedure must be reconsidered since the recursion coefficients 
Qn become double- valued in the large-n and N limit. A working hypothesis is to assume 
the following asymptotic form^^ 



Qn = Qc 



+ Ar-i/3(_)-/(t) + Ar-2/3(^(t) + 0(Ar-i), t = Ar2/3 ("i - ^V (55) 



Then the function A{X) in the limit ( |MD is modified by a constant term involving /'(O) 
and 5'(0), which are the same order 0{N~~'^^^) as A^ under the limit (^4]), 



AiX) = q„,{-2ql + ql + ql^, + A^) ~ N-^'^ {2g{0) - (-)"/'(0) + z'). (56) 

The numerical values of /'(O) and ^'(0) are determined by solving Painleve II equation 
(the microscopic limit of the Heisenberg algebra [Q, P] = 1) under an appropriate 
boundary condition^"^ . We may also have g run toward (?c = 1/4 with 
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N^/^{g-g,): fixed (57) 



while retaining the macroscopic spectrum intact. In this case the functions f'{t) and 
g{t) ought to be evaluated at t = 7 instead of t = 0. These numerical values enter the 
resulting second order linear differential equation for the wave function tpn (the counter- 
part of (^^). Thus the macroscopic spectral density p(A) is not sufficient to determine 
the microscopic correlation, but the latter depends upon the way how to approach 
the critical point. Microscopic universality holds in a weak sense for these cases, for 
Painleve II equation is universally derived from a class of the critical potentials. Details 
of the analysis will appear elsewhere^* . 
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